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r Integrals of the form 

nJ{Sb— cos X— cos y — cos z)dxdydz (1) 
\ Jo 

occur in problems that involve averages over the 
characteristic values, X, of the difference equation 



^^7l,m,n = ^yi 



(2) 



(where A^ is the Laplacian second difference oper- 
ator) with periodic boundary conditions 



yi,m,n y l-\-N ,171,71) yi,m,7i y 1 ,171 -{- N ,n f yi,7n,n~ 



-yi,7n,n + N' 

(3) 



Such averages arise in the theory of ferromagnetism, 
cooperative phenomena in solids, and critical fluctu- 
ations in liquids. 

In this paper, we deal with the tabulation of 



dxdydz 



(cos ic+cos ^+cos z) 



(4) 



[lyioT 6-^=/x=.oi(.oi)i. 

To obtain I{h) for values of /a from .01 to .8, we 
represent (4) as a power series in /x by the use of the 
following relations [1] : 

n(cos x+cos 7/+C0S zY dxdydz 
' JO 

=[^ir''""*r]...=''[£i«''>r]... 

where /o(<x) is the Bessel function of order zero for 
imaginary argument, which is given by the series 



XiT 



^"('^^-g-cio^ 



(6) 
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We obtain 



/(6)=f g(i)'"(2m)!c2„M'- (7) 

where the C2m are defined by 

[h{x)f = j: C,mX'^) (8) 



a tabulation of C2m (m = 0, 1, . . . , 20) is given in 
table 1. 

Table 1.« Coefficients of \h(x)f = y^.C2mX'^'" 



Numbers in parentheses stand for the number of zeros between the decimal 
point and the first significant figure. 



2m 


Com 





1. 




2 


. 75 




4 


. 23437 5 




6 


. 04036 4583 




8 


. 00433 34960 


9375 


10 


.00031 55517 


57812 5 


12 


. (4)16562 85038 


70081 02 


14 


. (6)65556 99977 


74500 42 


16 


. (7)20252 65474 


70871 28 


18 


. (9)50173 89455 


55226 63 


20 


. (10)10187 54486 


63855 90 


22 


. (12)17259 66354 


77837 00 


24 


. (14)24766 51388 


74193 16 


26 


. (16)30484 43411 


27630 85 


28 


. (18)32538 95366 


70253 00 


30 


. (20)30405 14661 


05592 65 


32 


. (22)25078 69824 


60377 59 


34 


. (24)18393 02011 


96660 07 


36 


. (26)12073 02467 


25867 06 


38 


. (29)71338 33750 


38164 30 


40 


. (31)38146 15500 


21396 56 



^Although the number of significant figures given here is many more than was 
required in the calculations for table 2, there may be some general interest in 
having these coefficients of [h(x)]^ available to many places. 

At most the first 21 terms of (7) were used to 
obtain 80 values of I{b) correct to five decimal places. 
The convergence of (7) is too slow to be used as a 
practical means of obtaining values of 1(b) beyond 

The remaining values of 7(6) given in the table 2 
were obtained by numerical integration. Since the 
integrand of (4) is infinite for x=y=z=0 and 6 = 1, 
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Table 2. 7(6) 




fl=b-' 


1(h) 


n = b-' 


Kb) 


.01 


0. 00333 


.51 


0. 17831 


.02 


. 00667 


. 52 


. 18219 


.03 


. 01000 


.53 


. 18610 


.04 


. 01334 


.54 


. 19003 


.05 


.01667 


. 55 


. 19400 


.06 


. 02001 


.56 


. 19799 


.07 


. 02335 


.57 


. 20202 


.08 


. 02670 


.58 


. 20607 


.09 


. 03004 


.59 


.21017 


. 10 


. 03339 


.60 


. 21429 


. 11 


. 03674 


.61 


. 21846 


. 12 


. 04010 


. 62 


. 22266 


. 13 


. 04346 


.63 


. 22691 


. 14 


. 04682 


.64 


. 23120 


. 15 


. 05019 


.65 


. 23553 


. 16 


. 05356 


. 66 


. 23991 


. 17 


. 05694 


.67 


. 24434 


. 18 


. 06033 


.68 


. 24882 


. 19 


. 06372 


.69 


. 25335 


. 20 


.06712 


.70 


. 25794 


. 21 


. 07052 


.71 


. 26259 


. 22 


. 07394 


. 72 


. 26730 


. 23 


. 07736 


.73 


. 27208 


. 24 


. 08079 


.74 


. 27694 


.25 


. 08422 


.75 


. 28186 


.26 


. 08767 


.76 


. 28687 


. 27 


.09113 


.77 


. 29196 


.28 


. 09459 


.78 


. 29714 


.29 


. 09807 


. 79 


. 30242 


. 30 


. 10156 


.80 


. 30781 


.31 


. 10506 


.81 


. 31332 


. 32 


. 10857 


.82 


. 31894 


.33 


. 11209 


.83 


. 32469 


.34 


. 11563 


.84 


. 33057 


. 35 


. 11918 


.85 


. 33662 


. 36 


. 12274 


.86 


. 34284 


.37 


. 12632 


.87 


. 34925 


.38 


. 12992 


.88 


. 35589 


.39 


. 13352 


.89 


. 36277 


.40 


. 13715 


.90 


. 36993 


.41 


. 14079 


.91 


. 37741 


.42 


. 14445 


. 92 


. 38526 


.43 


. 14813 


.93 


. 39355 


.44 


. 15183 


.94 


. 40238 


.45 


. 15555 


. 95 


.41189 


.46 


. 15929 


.96 


. 42229 


.47 


. 16304 


.97 


. 43395 


. 48 


. 16683 


.98 


. 44758 


.49 


. 17063 


. 99 


. 46503 


. 50 


. 17446 


1.00 


. 50546 



some modification was necessary before numerical 
integration could be applied. We make use of (4) 
in the form : 



7(6)=1 r^^-^ K (-j-^ )dx, 

ttVo 36 — cos X \36— cosx/ 



(9) 



where K is the complete elliptic integral of the first 
kind defined by 






dd 



^/l—u'sin^^ 



(10) 



In carrying out the integration, the integrand in 
(9) was expressed as the sum of two parts. The 
first part was the integrand diminished by an expres- 
sion that removed the singularity and made the 
numerical integration feasible. This expression was 
so chosen that it could be integrated in closed form. 

As an alternative to the foregoing method, the 
following procedure was employed as a check on 
some of the results. If we substitute the series 
expansion for K{u) in (9) and invert the order of 
summation and integration, we obtain ultimately 



7(6)=(95^-l)-i+Sri|: 



n = l L 




where ^n = [d"/^^(36)] {[9b'-l]-'^}. We note that g^ 
satisfies the following recurrence relation 

{9b'-l)gn+i + Sb(2n+l)gn+n'gn-i=0. (12) 
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